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During the past century, there has been considerable discussion and analysis of the motion of a 
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point charge in an external electromagnetic field in special relativity, taking into account "self-force" 
effects due to the particle's own electromagnetic field. We analyze the issue of "particle motion" in 
classical electromagnetism in a rigorous and systematic way by considering a one-parameter family 
of solutions to the coupled Maxwell and matter equations corresponding to having a body whose 
(3JT)! charge-current density J a (X) and stress-energy tensor T a i ) (X) scale to zero size in an asymptotically 



self-similar manner about a worldline 7 as A — > 0. In this limit, the charge, q, and total mass, m, of 
the body go to zero, and q/m goes to a well defined limit. The Maxwell field F a i,(\) is assumed to 
be the retarded solution associated with J a (A) plus a homogeneous solution (the "external field") 
that varies smoothly with A. We prove that the worldline 7 must be a solution to the Lorentz 
force equations of motion in the external field F a b(X = 0). We then obtain self-force, dipole forces, 
and spin force as first order perturbative corrections to the center of mass motion of the body. We 
believe that this is the first rigorous derivation of the complete first order correction to Lorentz force 
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motion. We also address the issue of obtaining a self-consistent perturbative equation of motion 
associated with our perturbative result, and argue that the self-force equations of motion that 
have previously been written down in conjunction with the "reduction of order" procedure should 
provide accurate equations of motion for a sufficiently small charged body with negligible dipole 
moments and spin. (There is no corresponding justification for the non-reduced-order equations.) 
We restrict consideration in this paper to classical electrodynamics in flat spacetime, but there 
should be no difficulty in extending our results to the motion of a charged body in an arbitrary 
globally hyperbolic curved spacetime. 
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I. INTRODUCTION 



In classical electrodynamics in special relativity, Maxwell's equations 

VF P = 4tt J M (1) 
V^ P ] = (2) 

are most often considered in the context where the charge current source J u is a specified 
function of the spacetime coordinates = (t,x % ). Provided only that J v is conserved, 

v v r = o , (3) 



Maxwell's equations have a well posed initial value formulation. Furthermore, since 
Maxwell's equations are linear in F^ u , it makes perfectly good mathematical sense to al- 
low J v be a distribution and to seek corresponding distributional solutions for F^ v . In 
particular, it makes perfectly good mathematical sense to consider point particle sources, 
i.e., distributional charge- current sources of the form (in global inertial coordinates 1 ) 

r = qu v ^\ x i -z\t)Y^ (4) 

where m m is the unit tangent (i.e., 4- velocity) of the worldline defined by x l (t) = z l (t), and r 
is the proper time along this worldline. Solutions to Maxwell's equations with an arbitrary 
point particle source of this form exist provided only that this worldline is timelike and that 
q is constant along the worldline (as required by conservation of J u , eq.([3])). 

Another problem that is commonly considered in classical electrodynamics is the motion 
of a point particle in a given external field. Here one prescribes a Maxwell field F^a as a 
function of spacetime coordinates and one postulates the Lorentz force equation of motion 
for the point charge, 

ma a = qF^fu , (5) 

where a a = w^V @u a is the 4-acceleration of the worldline. The Lorentz force equation is 
simply a system of second order ordinary differential equations for determining the worldline 

1 In arbitrary coordinates the right side would contain an additional factor of (— g)^ 1 ^ 2 (see eq. (l29|) below). 
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to which u a is the unit tangent. It is a well posed equation that admits a unique solution for 
the motion of the particle for any specified F^g (which need not be a solution to Maxwell's 
equations) and any specified initial conditions of the particle. 

Since Maxwell's equations are well posed for any specified point particle source and point 
particle motion is well posed in any specified Maxwell field, one might expect that the 
self-consistent coupled system, eqs.flT])-© and (jlj), together with 

ma a = qF a/3 u f3 , (6) 

would also be well posed. (Here, eq. flSj) differs from eq. (jSj) in that has been replaced by 
the full Maxwell field F a p, which includes the contributions from the "self- field" of the point 
charge.) However, it is easily seen that this is not the case: The full Maxwell field F a p is 
necessarily singular on the worldline of the particle, and, consequently, eq.® is ill-defined. 
Equations (jT])- (EJ) and (jlj) together with eq.([6j) simply do not make mathematical sense. 

Nevertheless, some authors — most notably Dirac |l[ — have attempted to proceed based on 
formal versions of these equations. However, in addition to being ad hoc, these treatments 
necessarily encounter mathematical and physical difficulties, such as the famous "infinite 
mass renormalization" . In this picture, the infinite electromagnetic self-energy of a point 
charge is compensated by a negatively infinite material mass, rendering the combined (ob- 
servable) mass finite. We see no justification for postulating that material mass could be 
negative (or infinite). Rather, as long as matter satisfies positive energy conditions, it is 
clear that a true point charge really would have an infinite total (observable) mass, and is 
therefore a physically unacceptable object. In short, the attempt to couple classical matter 
to electromagnetism via the point particle hypothesis (jlj) and the Lorentz force law (|3j) fails 
on both mathematical and physical grounds 2 . 

Thus, it seems clear that the point particle charge-current (jlj) and Lorentz force law 
(J5j) cannot be fundamental. Rather, matter must be modeled as a continuum, with these 
equations (as well corrections, including self-force effects) emerging as an approximate de- 
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the motion of "small bodies". In fact, much of the earliest work on the self-force 



0, Q, 13 



effectively followed this approach. However, a fundamental difficulty arises 



2 Similarly, in general relativity, Einstein's equation does not make mathematical sense for point mass 
sources [2j. Physically, matter would collapse to a black hole before a point mass limit could be achieved. 
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in this approach. If one takes a limit to zero size in the straightforward way (i.e., at fixed 
charge and mass), one will simply recover the point particle description and its associated 
difficulties, such as infinite self-energy. But if one keeps the body at finite size, then all of 
the body's internal degrees of freedom will affect its motion, and the description of motion 
will depend in a complicated way on the details of the structure and initial state of the body. 

This latter difficulty is commonly dealt with by postulating a rigid body, thereby elim- 
inating the internal degrees of freedom. However, the notion of the "rigidity" of a body is 
incompatible with special relativity unless the body follows the orbits of a Killing field. Thus, 
rigidity places extremely restrictive a priori constraints on motion in special relativity and 
is an unacceptable starting point for a general analysis of the motion of a body. A proper 
continuum model must correctly allow for relativistic deformations. Mathematically, one 
should consider only models such that the self-consistent, coupled, Maxwell-charged- matter 
equations admit a well-posed initial value formulation, such as a charged fluid or a charged 
elastic solid (see |6|] for further discussion). 

However, our interest here is in deriving "universal" properties of the motion of small 
bodies that do not depend upon the details of any particular matter model. Thus, we wish 
to assume only the existence of a matter stress-energy tensor T*^, which couples to the 
electromagnetic field via conservation of total stress-energy, 

V (T£ + T™) = , (7) 

where the electromagnetic stress-energy tensor is given by 

^ ( W - \g,uF aP F^ . (8) 
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We will refer to the coupled system of eqs.([I]|3]) and (171) as the Maxwell and matter equations. 
All realistic bodies composed of continuum matter and charge-currents should satisfy the 
Maxwell and matter equations, and our results will be derived solely from these equations. 3 
Note that the Maxwell and matter equations are nonlinear in the electromagnetic field. 



3 These equations were also used in [7( as the starting point in the self-force context. Very recently, this 
work was improved Q in a non-perturbative framework that provides a complete formalism for describing 
the motion of extended bodies, including self-field effects. 
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Since the motion of a finite-size body will depend on its detailed composition as well as on 
the details of its internal states of motion, we must consider a limit where the body is shrunk 
down to zero size in order to have a chance at obtaining a simple, universal description of 
its motion. As mentioned above, normally such a limit is taken at fixed charge and mass, 
thereby reproducing the same difficulties as arise if one attempts to work directly with 
point particles. The main new idea in the present work is to consider a mathematically 
precise limit — based on suitable one-parameter families of exact solutions to the Maxwell 
and matter equations — that will enable us to avoid these difficulties and thereby to derive 
self-force effects in a mathematically rigorous manner. Following the basic ideas of |s| (which 
were formulated in the context of general relativity), we consider a modified point particle 
limit, wherein not only the size of the body goes to zero, but its charge and mass also go to 
zero. More precisely, we will consider a limit where, asymptotically, only the overall scale 
of the body changes, so, in particular, all quantities scale by their naive dimension. In this 
limit, the body itself completely "disappears", and its electromagnetic self-energy goes to 
zero. We will show in section II II I that the worldline, 7, to which the body shrinks down 
must be a solution to the Lorentz force equations of motion. (As we shall also see in section 
IHIt the electromagnetic self-energy of the body makes a finite, non-zero contribution to the 
mass, m, that appears in the Lorentz force equation.) Self- force effects (as well as dipole 
force and spin force effects) then arise as a first order perturbative correction to the center of 
mass motion of the body. We emphasize that our approach removes from any fundamental 
status both the point charge source (jlj) and the Lorentz force law (jSJ). Rather, the notion of 
a "point charge" and the Lorentz force equation that it satisfies will be derived as a lowest 
order description of the electromagnetic field and bulk motion of a small (but extended) 
body 4 whose exact evolution is governed by the Maxwell and matter equations (JTH3J), and 

Our derivation will enable us greatly clarify what is perhaps the most famous difficulty 
of the electromagnetic self-force problem: "runaway solutions". The output of the classic 

n n n 

Abraham [3[, Lorentz [4J and Dirac |l[ analyses (as well as of many other analyses) is the 



4 Of course, the body of most interest to experimental phenomena, the electron, is an intrinsically quantum- 
mechanical object. Our results will apply to electrons only to the extent that a classical description of 
the motion of an electron can be justified. 
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following equation: 



2 

ma a = qF°fv,P + -q 2 {g a p + u a u p ) a , (9) 



usually known as the Abraham-Lorentz-Dirac (ALD) equation. The second term on the 
right side of this equation corresponds to the "self-force" on the charge resulting from its 
own electromagnetic field. This term makes the ALD equation third-order in time, rendering 
the status of its initial value formulation — and consequently its role in predicting actual 
particle motion — at best unclear. Furthermore, eq. ([9]) has solutions whose acceleration 
grows exponentially in time — behavior that is not observed in nature. Clearly, eq. ([9]) is 
disastrously pathological as it stands. 

Historically, there have been three basic responses to these problems. The first response 
views eq. (jHJ) to be a genuine prediction of classical particle electrodynamics, and seeks a 
resolution of the runaway difficulties by a treatment of the problem within the framework of 
quantum mechanics or quantum electrodynamics [10(] . If so, then classical electrodynamics 



would be seriously deficient. The second response [1[ holds that eq. J 9]) may be salvaged by 
introducing a set of rules designed to eliminate runaway behavior However, such rules 
are necessarily acausal in nature and would appear to be incompatible with the underlying 
causal behavior of classical electrodynamics. The third response maintains that the equation 



should be modified by the "reduction of order" procedure [111]. This procedure is based on 
the observation that eq. ([9]) is only expected to be valid to order q 2 . Therefore, if we replace 
the or term in eq. (Q by the value it would take for Lorentz force motion in the background 
field F^p, we should get an equation that is "equally valid" to order q 2 . Thus, to order q 2 , 
we may replace a a on the right side of eq. (Q by 

= ^K>" (io) 



m 



and we may replace a a by 



a fi = u% (2-Fjgu° 



U <VV 7 (±Fp) + (±Fff) ^-F^ s u 5 (11) 
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The resulting equation 



ma r 



m 



m y / m 
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is known as the "reduced order" Abraham-Lorentz equation. This equation is free from 
the "runaway" solutions that plague eq.([9]), and has a standard, second-order initial value 
formulation. Thus, eq. (fl2l) appears to be satisfactory on both mathematical and physical 
grounds. However, the reduction of order procedure has the appearance of being rather ad 
hoc, and if eq. is, in some sense, more fundamentally correct, it is difficult to justify 
replacing this equation with eq. (fT2l other than on purely pragmatic grounds. 

Our analysis sheds considerable light on this issue. We will prove in section [IV] that the 
first order deviation from Lorentz force motion is given by 5 

5[ma a ]=5[qFlfuP] + {g£ + u a uP) \l q ^R- ap - l -Q^V p Ff (a?S lP - 2u 5 Q\ p F$) }. 

(13) 

In this equation 8[. . .] refers to the first perturbative correction of ... , whereas all other 
quantities are lowest non-vanishing order. The parameters q, m, Q a/3 , and S a @ are, respec- 
tively, the charge, mass, 6 electromagnetic dipole moment, and spin of the body. Here the 
electromagnetic dipole Q a @ can be decomposed in terms of the electric dipole moment p a 
and magnetic dipole moment fi a by 



(14) 



with u a p a = u a fi a = 0, and the spin tensor S a/3 can be expressed in terms of the spin vector 
S a by 

S af3 = e *PlS u ^ Ss _ ( 15 ) 



5 We also obtain equations for the time evolution of S a b and dm (see eqs. (|135|) and (|136[) below). However, 
the electromagnetic dipole Q ab does not have any associated conservation law and can be "changed at 
will by the body itself" , so there exists no universal evolution equation for this quantity. Similar remarks 
should hold for all higher moments (of mass and of charge), although these moments would show up only 
at higher orders in our perturbation expansion. 

6 The perturbed mass, 8m, includes a contribution from the interaction energy of the electric dipole with 
the external field (see eq. (|131[) below). 



7 



The first term on the right side of equation (|T3j) is the first order correction to the ordinary 
Lorentz force law. The first term in braces corresponds to the ALD self-force. The second 
term in braces is a relativistic form of the usual dipole forces familiar from electrostatics 
and magnetostatics. The remaining terms are additional dipole and spin forces relevant in 
nonstationary situations. A non-relativistic version of these force terms is given in eq. (|143j) 
below. 

In the case of a body with negligible electric and magnetic dipole moments and negligible 
spin, only the first two terms of equation (fT3|) remain. In this case, we see that the form of 
the ALD equation ([9]) is recovered as the leading order perturbative correction to Lorentz 
force motion for a sufficiently small charged body. However, since the acceleration a M present 
in the ALD term refers to the "background motion" (i.e., a solution to the Lorentz force 
equations of motion in the background electromagnetic field), no runaway behavior occurs 
in the context of perturbation theory. 

Although ( fT3l) has been rigorously derived as the leading order perturbative correction 
to Lorentz force motion, small deviations from Lorentz force motion can accumulate over 
time. Eventually, the deviation from the zeroth order motion should become large, so one 
would not expect that solutions to (fT3l) would provide an accurate description of motion at 
late times. However, as argued in 9|, if the deviations from Lorentz force motion are locally 
small, it may be possible to find a "self-consistent perturbative equation" based upon eqs.flS]) 
and (fT3l) that provides a good, global-in-time description of the motion. We shall argue in 
section|V]that in the case of a body with negligible electric and magnetic dipole moments and 
negligible spin, the reduced order Abraham-Lorentz equation ( 1T21) provides such a desired 
self-consistent perturbative equation. This provides a solid justification for using eq. (fT2]) to 
describe the motion of sufficiently small charged bodies (with negligible dipole moments and 
spin). There is no corresponding justification for using the ALD equation (J9]) to describe 
the motion of charged bodies. 

Our notation and conventions are as follows. Greek indices fi, u, .. will refer to spacetime 
coordinate components, whereas mid-alphabet Latin indices will refer to spatial co- 

ordinates only. A "0" will denote the time component. Early-alphabet Latin indices a, b, ... 
will refer to abstract spacetime indices 12|. (Typically, a,b, ... will be used in diffeomor- 
phism covariant expressions, whereas fi,u, ... will be used for expressions that take a special 
form in the coordinates being used.) Our (flat) metric g a b has signature (—1,1,1,1). We 
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normalize the Levi-Civita symbol e a b C d so that eoi23 = +1- 

II. ASSUMPTIONS 

We consider a one-parameter family of charged bodies described by a charge-current 
vector J M (A) and a matter stress-energy tensor Tj£(\). The electromagnetic field -Fj^(A) 
is assumed to satisfy Maxwell's equations (jTJ and (j2J), the charge-current is assumed to be 
conserved, eq.(E]), and the total stress energy tensor is assumed to be conserved, eq.flZj). Since 
we seek relations that hold universally for all bodies, we specify no additional "constitutive" 
relations, and impose only eqs.([T])-([3]) and ([7]) in our analysis. 

As discussed in the previous section, we wish to consider a limit where the body shrinks 
down to zero size, but in such a way that the mathematical (and physical) inconsistencies 
associated with point particles are avoided. Following the basic ideas of [9], this can be 
done by considering a limit where not only the size of the body goes to zero, but its charge 
and mass also go to zero in such a way that, asymptotically, the body continues to "look 
the same" except for its overall scale. In [9(, this idea was implemented in a relatively 
indirect way by postulating the existence of a "scaled limit" of the exterior field of the body 
together with the existence of an "ordinary limit" and a "uniformity condition" on these 
limits. In the context of general relativity, it was necessary in [9] to formulate the conditions 
on the one-parameter family in terms of the exterior field because we wished to consider 
strong field objects (such as black holes), which do not admit a simple description in terms 
of a stress-energy source. However, in the present context of classical electrodynamics in 
Minkowski spacetime, all of the relevant information about the body and the electromagnetic 
field is contained in the specification of the charge current density J M , the matter stress- 
energy tensor Tr£, and "external" electromagnetic field, i.e., the homogeneous solution of 
Maxwell's equations that describes the difference between the actual electromagnetic field 
and the retarded solution associated with J M . Consequently, we will formulate our conditions 
in terms of these quantities. 

We wish to consider a one-parameter family {F fJjU (X), J M (A), T^(X)} of solutions to 
eqs. (p])- (ED and (0) having the property that the worldtube, W(A), containing the supports 
of J M (A) and T*£(X) shrinks down to a timelike worldline 7 as A —>■ 0. Furthermore, we wish 
to impose the condition that this "shrinking down" asymptotically corresponds merely to a 
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"change of scale". To get a feeling for what this entails, fix a time, to, on 7, and consider 
the charge density, p(X,x l ), on the hyperplane orthogonal to 7 at to given by 

p(A,x*)=A- 2 p(^y^) (16) 

where p is a smooth function on R 3 that vanishes for r > R for some R > 0. Then p(A, x % ) 
vanishes for r > XR, i.e., the support of p shrinks to z l as A — > 0. It also is clear that for the 
one-parameter family, eq.f fTBl) . the charge density "retains its shape" exactly as it shrinks 
to the origin. Finally, the overall scaling of A -2 in eq. (116ft ensures that the total charge also 
goes to zero proportionally to A. Thus, eq. ffTB"]) describes the kind of "scaling to zero size" 
that we seek. 

However, since a body undergoing non-uniform motion cannot remain "rigid", we must 
allow for time dependence in the manner in which the charge-current and body stress-energy 
scale to zero size, particularly in view of the fact that the body may find itself in a different 
external field at different times. In addition, the exact scaling represented by eq.f fTB"]) is too 
strong an assumption; in particular, it could hold, at most, for a particular choice of time 
slicing and a particular choice of spatial coordinates on these slices. Furthermore, although 
we could always impose eq. ffTB"]) exactly as an initial condition at one time, there is no reason 
to believe that dynamical evolution would preserve this exact scaling at later times. Instead, 
it would be much more reasonable to require that the type of scaling represented by eq.ffTBl) 
hold only asymptotically as A — > 0. 

In order to formulate our conditions on J M (A) and T/J5(A), it is useful to adopt a spe- 



cific choice of coordinate system, namely Fermi normal coordinates (see, e.g., 13]). These 
coordinates (T, X 1 ) are defined with reference to a timelike worldline, 7. The surfaces of 
constant T are orthogonal to 7, and the X 1 measure spatial distance within those surfaces. 
With this coordinate choice we may state our conditions on J M (A) and T^(X) as follows: 

• (i) We require that there be a smooth, timelike worldline 7 such that in Fermi nor- 
mal coordinates, (T, X 1 ) , based on 7, we have J^(A, T, X*) = A~ 2 J M (A,T, X i /X) and 
T^(A, T, X*) = A- 2 f*J(A,T,X7A), where f M and J are smooth functions of their 
arguments. Furthermore, the supports of T M and J in the last variables — i.e., the 
slots into which X l /X have been inserted — are contained within a radius R, so that 



the J M (A) and Tf£(\) are nonvanishing only within a worldtube, W(A), of radius XR 

10 



about 7. 

Note that the smooth dependence on both T and A permitted for J M and T^ v in this condition 
addresses the concerns of the previous paragraph. 

We now present these assumptions in general coordinates. Let (t, x l ) denote any coor- 
dinates that are smoothly related to Fermi normal coordinates in the neighborhood of 7 in 
which the Fermi normal coordinates are defined and are such that V M t is timelike. In these 
new coordinates, the worldline 7 will be given by an equation of the form x % = z l (t). Since 
the Fermi coordinate X 1 is a smooth function that vanishes on 7, we must have 

where hij is smooth. Writing T = T(t,x ), writing 

A 

and substituting into the formula for the charge current density in the new coordinates, we 
find that the components of the charge current density in the new coordinates take the form 

J"(A, t, x') = A~ 2 J /M (A, *, [x* - z l {t)}/\) , (19) 

where J' M is a smooth function of its arguments. Furthermore, since V M t is timelike, each 
surface of constant t must intersect W(A) in a compact set, which implies that J /M is of 
compact support in the last variables. Note that, in particular, the form eq. (lT9~]) holds in 
global inertial coordinates. By reversing the argument, it is easily seen that assumption 
(i) on J M (A) as stated in Fermi normal coordinates is equivalent to eq. (|T9|) holding in some 
coordinates with V M t timelike, together with the requirement that the support of J M (A) be 
contained in a worldtube of proper distance XR from 7. Similar results hold, of course, for 
T^(\,t,x>). 

Our assumption about the electromagnetic field F^ U (X) is as follows: 
• (ii) We have = Fj%? + F^ lf , where F^ lf is the retarded 7 solution of Maxwell's 

7 It may appear that we are introducing a time asymmetry in our assumptions by choosing F^} { to be the 
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(17) 



(18) 



equations with source J M (A) and Fj^ is a homogeneous solution of Maxwell's equation 
that is jointly smooth function of A and the spacetime point. 

Assumptions (i)-(ii) about our one-parameter family together with eqs. ([!])- (J3]) and (J7J) 
constitute the entirety of what we shall assume in this paper. 

In Appendix |Aj we show that in arbitrary smooth coordinates {t,x l ), the retarded solu- 
tion, F^ lf (A), with source of the form eq.(lT9l takes the form 

F;f(\,t,x*) = X-'F^t^x* - z\t)]/\) , (20) 

where F is a smooth function of its arguments. 

Finally, we relate the assumptions about our one-parameter family made here to the 
assumptions made in 9]. As previously mentioned, the assumptions of [9] were formulated 
entirely in terms of the behavior of the "exterior field" of the body, i.e., the behavior of 
the spacetime metric outside of the body. Assumption (i) of [9j ( "existence of the ordinary 
limit") corresponds here to the requirement that F^ be jointly smooth in (X,t,x t ) except 
on 7 and that F^ U (X = 0) be smoothly extendible to 7. That this condition holds follows 
immediately from our assumption (ii) together with the fact that goes to zero at fixed 
(i, x l ) with x % 7^ z' L {t) as A — * 0. The analog of assumption (ii) of [9] ("existence of the 
scaled limit") follows immediately from assumption (ii) and eq.( l20|) . We will consider this 
scaled limit in detail in section IIVI below. Finally, the analog of the "uniformity condition" 
(iii) of 9J is as follows: Define a = r = yfYl^x 1 ~~ z l (t)] 2 and define f3 = X/r. Write F^ v 
as a function of (a, (3, t, 9, 0). The desired condition is that at fixed (t, 9, 0), XF^ U is jointly 
smooth in (a, (3) at (0,0). Since XF^ is easily seen to satisfy this condition, the analog of 
the uniformity condition of [9J will hold here provided that F^ u is smooth in (a, (3) at (0, 0) 
at fixed (t,9,(p). In appendix [A| we show that not only is F^ v smooth in (a,/3), but is of 
the form of (3 2 times a smooth function of (a,/3). 

Thus, the assumptions we have made above imply analogs of the assumptions made in 
9|. Since our assumptions here are quite simple and straightforward as compared with the 
assumptions made in j^, this may be viewed as providing justification for the assumptions 



retarded solution. However, choosing F^} 1 to be the advanced solution results in an equivalent set of 
assumptions, since the retarded minus advanced solution is a homogeneous solution that can be shown to 
be smooth in (X,t,x l ). 
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made in |9j. Since the assumptions (i)-(iii) of 9] are essentially what is needed to justify 
the use of "matched asymptotic expansions" [13l. Il4|. we have thus effectively also provided 
justification for the use matched asymptotic expansions. 



III. THE FAR-ZONE LIMIT AND LORENTZ FORCE MOTION 

As in [91], as A — > it will be useful to consider both a "far- zone" limit (wherein x^ is 
held fixed) and a "near zone" limit (wherein x^ = ([£ — to]/\, [x % — z l (to)]/\) is held fixed). 
Roughly speaking, the "far zone" limit corresponds to how the body and the electromagnetic 
field appear to an observer at a fixed radius from 7 as the body shrinks down to 7. The near 
zone limit describes the appearance to an observer who "follows the body in" toward 7 at 
time to — and correspondingly rescales units — as the body shrinks. We use the terminology 
"far zone" and "near zone" because of the close correspondence to the way these terms 
are used in matched asymptotic expansion analyses. However, we will not need to do any 
"matching" in our analysis, and any of our calculations can be done in either the "far" or 
"near" zone pictures. In this section, we treat the far- zone limit. We will show that at first 
order in A, the far zone description of the body is that of a point particle, and that the 
worldline 7 satisfies the Lorentz force equation of motion ((Sj). Interestingly, we will also see 
that the electromagnetic self-energy of the body makes a non-zero, finite contribution to the 
particle's mass. 

It is clear from assumption (i) of the previous section that J' 1 (A) and T^(X) as well as 
all of their A-derivatives go to zero pointwise as A — > and any fixed x^ ^ 7. On the other 
hand, J^(A) and T*J(A) and their A-derivatives do not approach any limit at all for x M G 7. 
Nevertheless, J M (A), Tf^{X) and their A-derivatives have well defined limits as distributions 
as A — > 0. To see this, we recall that associated to any locally L l function H(t, x l ) on 
spacetime is the distribution, Dh, defined by the following action on a smooth test function 
f(t,x l ) of compact support, 




(21) 
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Thus, viewed as a distribution, the charge- current J^{X) is given by 



Dj(\)[f]= I r{\t,x % )f^x 1 )^ dtd 3 x 

A~ 2 >(A,t, [x* - z i {t)]/X)f IJL {t,x i )y/=g~ dtd 3 x 
= Xj J tl {X,t,x i )f ll {t,z i {t) + Xx^y/^g dtd 3 x 
= xj dt^Mt,z\t)) J (PxJtoJffrt,?) + 0(A 2 ) . (22) 

(Here we work in arbitrary smooth coordinates (t,x i ) and have used eq. (|T9|) in the second 
line; we also have dropped the prime on J M .) Thus, we see that, viewing J M (A) as a one- 
parameter family of distributions, we have 

J (0)M = lirnJ^A) = (23) 

jd), = lim JL j^ X ) = J»(t) S{xl (24) 
a^o dX v 1 w dt 



with 



J»(t) = J >(A = 0,t,x l )^~ 3 d 3 x, (25) 



Higher-order derivatives of J M with respect to A will similarly have the form of distributions 
with support on 7. For example, j( 2 )^ will contain a delta function, which will contribute 
a correction to the charge of the body, along with a derivative of delta function, which will 
provide the leading order contribution to the electromagnetic dipole moment. In general, 
the nth order charge-current j( n ^ will have the form of the sum of derivatives of the delta 
function up to order n, where the coefficient of the mth order delta function is the {n — m)th 
correction to the (n — l)st moment. 

We can obtain considerably more information about the form of J^ 1 ^, eq. (124"]) . by making 
use of conservation of J M , eq.Q. For any test function /, we have 

J V u r{X)f^-gd i x = , (26) 

for all A > 0. Integrating by parts, differentiating with respect to A, and taking the limit as 
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A — > 0, we obtain 



J J f "(t)[W,f](t,z\t))dr = 0, (27) 



for all test-functions /. Now, for any i — 1, 2, 3, choose / to be a test-function of the form 
/ = X l c(t)h(Xi), where X % is a Fermi normal spatial coordinate and c and h are of compact 
support with h(X^) = 1 in a neighborhood of the origin. Equation (12?]) then yields 

J J^t){V fl X l ]c(t)dr = 0, (28) 

which implies that the projection of j7 M (t) orthogonal to 7 vanishes for all t, i.e., J 11 is of 
the form J(t)u^. Using this result in equation ([2"T]) and integrating by parts, we find that 
j{t) = 0. We thus obtain 

jV» = qu ^-^-£M *L, (29) 

V^g dt 

where the (lowest-order) total charge q is given by 

q= / j°(\ = Q,t,x i ) y /(fy{Px (30) 



and is independent of time t. Thus, to first order in A, the charge current of the body is 
precisely the standard form of a structureless point charge moving on the worldline 7. 

The electromagnetic field F^ u is the sum of and F^ f . The external field Fj%? is 
smooth in A and satisfies the homogeneous Maxwell's equations at all A. Therefore, -^^(A) 
has a straightforward perturbation expansion in A, and all A-derivatives of Fj%? satisfy the 
homogeneous Maxwell's equations at each order in perturbation theory. 

From eq. (IA22I) of appendix we see that the self-field is of the form 

F »? =^(^,A/r,0,0) , (31) 
where, for any fixed (t, 9, 0), is smooth near in its second and third arguments. Using 
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this property, we may Taylor expand F^} { to any finite orders N, M, 



fJtU 



x N M / \ \ ' m 

F°?(t, r, M) = £ £ £ r n ~ r (^W*. 0, 0) + O (r" +1 ) + O ((A/r)^ +1 ) 

n=0 m=0 ^ ' 

AT M 

= £ £ A m+1 r"~ m - 2 (^) nm (t, 61, 4>) + O {r N+1 ) + O ((A/r) M+1 ) (32) 



n=0 m=0 



This gives the general form of the far-zone expansion of F^l near r = 



We can obtain explicit forms for the terms arising in the far-zone expansion of F 



self 
fj,u 



by using the fact that for all A, F^ v (A) is the retarded solution associated with the charge- 
current J M (A), so the n-order perturbative self-field Fj^'^ 1 ' is the retarded solution associated 
with the distributional source J^'^. In particular, the zeroth order self-field vanishes 

F self,(0) = (33) 

and the first order self-field F^t ' is just the standard Lienard-Wiechert solution associated 
with the point charge source, eq. (l24l) . In Fermi normal coordinates about 7, F^, lf ' takes 
the explicit form 

r seif,(i) [n% 1 / I , 1 A 

F i0 ^ = g j- + - \--( aj n')m - -Oij + 

+ ~(a,j'n?) 2 Tii + — (ajri')ai + -aja^rii + -a ni + -d, + 0(r) 

Ff>M = - q ±a [inj] + 0(r), (34) 

where n l = x l jr. 

We now turn our attention to the stress-energy tensor. By assumption (i), the matter 
stress-energy T^(A) has the same type of behavior as J M (A) as A — > 0, and thus it has an 
analogous multipole series. In particular, as a distribution, T^(\) — > as A — > 0, whereas 
to first order in A, the matter stress-energy tensor takes the form 

t m,(i) _ i- d T A/ m _ T M (+ \ b{x l - z l (t)) dr 

T ^ = 1™ d\ T ^ {X) ~ T ^ {t) 7=1 it ■ (35) 

The electromagnetic stress energy tensor is given by eq.flU), with F^ v = F^+F^ f . Thus, 
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the total electromagnetic stress-energy naturally breaks up into a sum of three terms, 

rpEM ry-iext I ry-icross , ry-iself 

where the notation should be self-explanatory. The stress-energy tensor Tfj^ of the external 
field behaves smoothly in (A, x^) and thus has a straightforward perturbative expan- 
sion in A. Since Ff^ is a solution to the homogeneous Maxwell's equations, we also have 
V /i T^ t (A) = for all A, so conservation of T^f holds at each order in perturbation theory. 

The behavior of the stress-energy tensor Tj^f of the self-field F^f as A — > is more 
delicate. Since F^f { '^ = 0, one might expect that the lowest order contribution to T*f { would 
be obtained by substituting XF^f^ 1 , eq. fjM]) . into eq.©. This would give an expression 
of the form A 2 times the usual stress-energy tensor of a point charge. This would suggest 
that the first order in A contribution to T^ lf vanishes, and the second order contribution is 
too singular on 7 to define a distribution. If this were the case, then our analysis would be 
plagued by the same type of infinite self-energy problems that occur in usual treatments of 
point particles. Fortunately, as we now shall show, a more careful treatment shows that to 
first order in A, T^ lf has a non-zero, well defined distributional limit. 

As previously mentioned in the Introduction, in Appendix |A] we prove that F^ lf is of the 
form (1201) . where F^ u is a smooth function of (X,t, [x l — z l (t)]/X). We also prove in Appendix 
[A] that if we define a = r = \fYl,[x l — z l (t)] 2 , define (3 = A/r, and write F^ u as a function 
of (a, P, t, 9, (f>), then, at any fixed (t, 9, (ft), F^ is of the form of f3 2 times a smooth function 
of (a, (3) at (0,0). Consequently, T^ lf is of the form 

T*f(X, t, x*) = A- 2 f lf(X, t, [x* - z\t)}/X) , (37) 

where is given in terms of F^ by eq. (jSJ) and is a smooth function of its arguments. Thus, 
(A) is of a form similar to that of J*(\) and Tff v {X) in our condition (i) of the previous 
section, except that it is not of compact support in x l . Nevertheless, from the properties of 
F^, we know that at any fixed (t,9,(p), T^ lf is of the form of /3 4 times a smooth function 
of (a, (3). If we view T^ lf (A) for A > as a distribution, the same manipulations as led to 
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eq. (!22|) now yield 

D Tse if[/] = A j T^ lf (A, t, x^f^it, z\t) + Xx^^dtd^x . (38) 

If T^ lf were of compact support in x l , we could straightforwardly take the limit as A — > 
inside the integral, as we did in eq. (l22l) . Nevertheless, using the facts that (1) the integrand 
is smooth in x\ (2) the test tensor field f^ u is of compact support in x % (and, thus, is of 
compact support in a) and (3) for a in a compact set, T^ lf is bounded by C/3 4 for some 
constant C, it is not difficult to show that for A < Aq 



\f;f(\,t,x l )\\r(t,z\t) + Ax 4 ) I < -j— (39) 



K 

H 



where K is a constant (i.e., independent of A and x l ). Since l/[f 4 + 1] is integrable with 
respect to y/—gd 3 x, the dominated convergence theorem then allows us to take the limit as 
A — > inside the integral in this case as well, and we obtain 

D TSC «[f} = \J dt^r(t,z l (t)) J d 3 x^f*f(0,t,x*) + O(\) , (40) 

i.e., as distributions we have 

^f(°) = limT- lf (A) = (41) 

A — >0 



T self,(l) = j. ^_ T sel{ (x) _ r self/ t \^ Z l (t)) <fr 



with 



= j f^(0,t,x l )^d 3 x . (43) 



Thus, despite the fact that the electromagnetic self-field is not confined to a world-tube at 
any A > 0, to first order in A, the electromagnetic self-stress-energy takes the form of a 
5-function on 7. However, it should be noted that, in contrast to the situation for J M and 
T^, at higher orders in A the support of T^ lf cannot be confined to 7. Indeed, although we 
have shown above that D TB cu[f], eq. (l38l) . is a C 1 function of A at A = 0, it does not appear 
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that it is smooth 8 in A at A = 0. 

There is nothing further that we can say beyond eqs. ()35l) and (j4"2l) about the individual 
forms of the material stress-energy and electromagnetic self-stress-energy to first order in A. 
However, defining the "matter plus electromagnetic self- field" stress-energy T^, 

T, u = + T*f, (44) 

we see from conservation of total stress-energy, (171), that for all A we have 

= -v% c f - v% c ™ ss 

= F^J\ (45) 

where we emphasize that T^, with no superscript, includes the matter and "electromagnetic 
self- field" contributions, as in equation (1441 . The same type of calculation that led from 
eq. fl26|) to eq. fl29|) then gives to first order in A that 

. 5(x l — z l (t)) dr , irxS 

Tg>(t) = mu,u v ^_ Ul — (46) 

where u M is the unit tangent to 7, and m is a constant (i.e., independent of t). In coordinates 
(t, x l ) such that the constant-t hyperplanes are orthogonal to u^, m is given by 

m = j (f™ (A = 0, t, x l ) + f s clf (A = 0, t, x 1 )) V9~3 d 3 x. (47) 

Furthermore, must satisfy 

mu'V^ = qu u F^{\ = 0, t, z\t)) . (48) 

We may summarize what we have just shown as follows: Consider any one-parameter 
family of bodies satisfying the assumptions of the previous section. Then, to first order in 
A, the description of the body is precisely that of a classical point charge/mass moving on 



A failure in obtaining a distributional series for T^° lf should not affect our analysis at higher orders, since 
T^f f enters our equations only in the form V'T^° lf . It is not difficult to show that V M T^ lf does have a 
distributional series with support on 7 of a form similar to that of J^ 1 and TM,. 
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a Lorentz force trajectory of the external field. Specifically, the charge- current of the body 
is given by eqs. (1291) and fl30|) . the sum of the body stress-energy and the electromagnetic 
self- field stress-energy is given by eq.lHol). and the world line is given by eq. lflgj) . It should be 
noted from eq. (HT|) that the electromagnetic self-energy contributes to the mass of the body, 
i.e., there is a finite "mass renormalization" of the body contributed by its electromagnetic 
field. It should be emphasized that the description of the body as a point charge/mass is 
an output of our calculations, that our derivation is mathematically rigorous, and that no 
infinities have arisen in the derivation of our formula (H7|) for the mass of the particle. 

IV. NEAR-ZONE LIMIT AND PERTURBED MOTION 

We now turn our attention to the first order in A corrections to the motion. In order to 
do so, we must first address the issue of what we mean by the "corrections to motion" . As 
A — > 0, the body shrinks down to a worldline 7, so at lowest order, the motion is described 
by 7. As shown in the previous section, 7 must satisfy the Lorentz force equation (1481) . But, 
at any A > 0, the body is of finite size, so in order to find the "correction" to 7 at finite A, 
we would need to have a "representative worldline" to describe the motion of the body. It 
would be natural to take this worldline to be the "center of mass". However, as we shall 
now explain, it is far from obvious as to how to define a "center of mass" . 

The difficulty arises from the fact that "electromagnetic self-energy" must be included 
in a definition of center of mass. Indeed, we have already seen in eq. (H7I) above that, at 
lowest nonvanishing order, electromagnetic self-energy contributes to the mass of the body. 
Since energy can be exchanged between the matter and electromagnetic self-field — only 
their sum, m, is conserved at leading order — it clearly would not make sense to omit the 
electromagnetic self-energy from the definition of the center of mass. However, if one defines 
the electromagnetic self-stress-energy T^f at finite A to be the electromagnetic stress-energy 
of the retarded solution then the stress-energy of radiation that was emitted by the 

body in the distant past would be included. Furthermore, integrals over the self-stress- 
energy defined this way — corresponding to that stress-energy's contributions to mass, center 
of mass, etc. — would diverge, on account of the slow, 1/r falloff of the retarded self-field. 
It is therefore clear on physical and mathematical grounds that the full retarded self-field 
may not be used to provide a definition of center of mass. On the other hand, it is far from 
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obvious as to how one can exclude the energy of the "emitted radiation" while including the 
energy of the "Coulomb self-field", since one cannot make a clean distinction between the 
two at finite times. 

n 

Recently, one of us [8j] has proposed an appropriate definition of the "center of mass" 
for charged bodies that properly takes into account the body's electromagnetic self-energy 
while not including radiation emitted in the distant past. It should be possible to use the 
"representative worldline" provided by this definition at small but finite A to define a notion 
of the perturbative corrections to 7 to all orders in A. However, in this paper, we shall 
be concerned only with the leading order correction to 7. As we shall see below, at this 
order, there is a very simple and straightforward notion of "center of mass" — including self- 
field contributions — that can be used to define the perturbed motion, and we will use this 
notion in our analysis. It can be shown that this notion agrees with the much more general 
definition of [8( to this order. 9 

Our analysis is most conveniently carried out in the "near zone" picture. Therefore, 
in subsection IIV Al we introduce the near-zone limit for the charge-current, stress-energy, 
electromagnetic field, and metric, and we determine the properties of the various near-zone 
fields. In subsection IIV we present our definitions of the body parameters. Finally, in 
subsection IIV CI we re-derive the background Lorentz force motion in the near-zone context, 
and then compute the perturbed motion. 

A. Near-zone Limit and Properties 

As analyzed in detail in the previous section, the ordinary ("far-zone") limit associated 
with our one-parameter family sees the body shrink down in size to an effective point particle 
description. Here we define a second, "near-zone" limit wherein the body remains at fixed 
size. To accomplish this, we measure spacetime distances (at A > 0) with a rescaled metric, 

g ab ee \- 2 g ab . (49) 

9 It was emphasized in Q that there is a class of reasonable definitions of generalized momenta and hence 
center of mass that one may adopt. The differences between the definitions begin at 0(X 2 ) in the pertur- 
bation series discussed in this paper, which is precisely the order at which our type of definition would 
break down. 
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The inverse metric g ab is related by g ab = \ 2 g ab . The Maxwell and matter equations, 
however, are not satisfied relative to the barred metric g ao . Therefore, we also introduce 
rescaled fields, 10 

J a = \ 3 J a (50) 

= T% (51) 
F ab = \- l F ah , (52) 

so that the Maxwell and matter equations, 

V b F ab = 4nJ a (53) 

V [a F % = (54) 

V a J a = (55) 

V b (^+^ M )=0, (56) 

are satisfied. Here the barred metric g ab is used to raise and lower indices, and V a is 
the derivative operator associated with the barred metric (which, however, agrees with the 
derivative operator associated with the unbarred metric). The electromagnetic stress-energy 
T™ is constructed from F ab in the usual way. 

To complete the description of the near-zone limit, we introduce at each time t the 
"scaled coordinates" (f, x 1 ) defined by 

t=*-=^, x ^-f (tQ) . (57) 

The limit as A — > at fixed (i, x l ) of a barred quantity in barred coordinates constitutes its 
near-zone limit. Physically, as A — > 0, an observer at fixed (i, x l ) approaches the spacetime 
point (t = t Q ,x % = z l (t )) along with the shrinking body, but since the observer uses the 
rescaled metric to measure distances, the space and time intervals between the observer 
and the point (t ,z l (to)) remain finite. Thus the near-zone limit has the interpretation 



10 Note that these rescalings are modified if indices are raised or lowered, since we use the barred metric 
to raise and lower indices for barred quantities and the unbarred metric to raise and lower indices for 
unbarred quantities. For example, we have J a = \J a . 
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of "zooming in" on the body (in space and time) at the same rate at which the body is 
shrinking. 

In order investigate the properties of the near-zone limit, it is convenient to specialize to 
coordinates where z l (t) = 0, and we shall do so for the remainder of this paper. In this case, 
near-zone coordinate components are related to far-zone coordinate components by simple 
expressions, 



g pp (X, t ; t, x*) = g^(X, t = t + At, x i = Ax*) (58) 

P(X, t ; i, x l ) = A 2 J M (A, t = t + At, x* = Ax*) (59) 

f^(A, t ; i, x l ) = A 2 T^(A, t = t + At, x* = Ax 1 ) (60) 

Fp P (X, t ; t, x*) = XF fMU (X, t = t + At, x* = Ax*) . (61) 



Note the presence of bars on the indices for the tensors on the left side of these equations 
(denoting components in barred coordinates) and the absence of bars on the indices of 
tensors on the right sides (denoting corresponding components in the corresponding unbarred 
coordinates). It then follows directly from our assumptions and eq. (120]) that 

P{X, t ; t, x*) = >(A, t + At, x*) (62) 
T£(A, t ; t, x l ) = T^(A, to + At, x l ) (63) 
Fp,p(X, t ; t, x l ) = F^(A, t + At, x l ) + XF°f(X, t + At, Xx l ) . (64) 

This gives the behavior of the near-zone fields in terms of the (smooth) "tilded" quantities 
of our assumptions. Equations (l62~!) - (l64l) make clear that in the near-zone limit the size of 
the body remains finite (as measured by the barred metric) and all properties of the body 
(as described by the barred fields) go to well-defined, finite limits. We also see that the 
effects of the external universe, as well as of the body's own future and past, reduce in 
importance and finally disappear at A = 0. This contrasts with the far- zone limit analyzed 
in the previous section, where the effects of the external universe go to well defined, finite 
limits, but the effects of the body reduce in importance and finally disappear at A = 0. 

Since the near-zone fields have now been seen explicitly to be smooth in A, i, x l , we 
know in particular that the near-zone perturbation series exists to all orders in A. We 
will denote the terms in this series in analogy with the notation used in the far-zone limit: 
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For a generic barred quantity /, we write f( n \x ,M ) = ^ lim A ^ J^f(\, x^). Thus it is our 
convention that the superscript (n) refers to A-derivatives at A = at fixed x M when the 
superscripted quantity is barred; otherwise, it refers to A-derivatives at fixed x M . Note that 
since the charge-current and material stress-energy are of compact support in x % at all A, 
this property also holds at each order in perturbation theory. 

It follows immediately from eqs. (l58]) - (l6i~]) that each near-zone quantity / depends on to 
and t only in the combination i + At, and that / must be a smooth function of to + At. 
This has two important consequences. First, the dependence of any nth order perturbative 
quantity p n ' on t can be at most nth order polynomial in t 

n 

f n \i,x l ) = J2 A % ) (x l )i m - (65) 

m=0 

In particular, all zeroth order near-zone quantities are stationary. Second, we have a simple 
relationship between t and to derivatives of p n \ 

dt J <9t 

Equations ( |65i) and ( 1661) are satisfied component-by-component by all the near-zone quanti- 
ties JP 1 , T~, F^ p , g^p. Equation fl66l) reflects the fact that the near-zone perturbation series 
at each time to are all defined with reference to the single ( "far-zone" ) one-parameter-family. 
Thus, perturbed near-zone quantities at different scaling times to cannot be specified inde- 
pendently but instead must be related by equation (ItjfiT) . We refer to this relation here, as 
in we did for the analogous condition in [9( , as a "consistency condition" . We will make use 
of the consistency condition in our calculations, below. 

We now turn to the relationship between near-zone and far- zone perturbative quantities. 
In the previous section, we saw that the far-zone expansion of J M was distributional in 
nature, with the (n + l)st order term, j( n+1 )^ in the expansion being described as a sum 
of multipoles (i.e., derivatives of 5-functions) up to order n. It can be seen that at t — t , 
for m < n, the mth order multipole (in the sense of an mth-derivative of a delta function) 
appearing in j( n+1 )^ contributes to the mth multipole moment (in the sense of a moment 
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integral) of J( n - m )A a t f = o. 11 Conversely, the multipole moments of the near- zone pth 
order charge- current J( p ^ correspond to moments appearing at pth and higher order in the 
far- zone expansion of J M . This is characteristic of the kind of "mixing of orders" of quantities 
appearing in far- zone and near-zone expansions. Similar results hold for the near-zone and 
far- zone expansions of T^. 

Consider, now, the near-zone external electromagnetic field 



»(A, t ; t, x l ) = AF c f (A, t + At, Ax*). (67) 



Since Ff^ is a smooth function of its arguments, we see that the perturbative expansion 
of .F^p* at order n can depend only polynomially on x l , with the degree of the polynomial 
being no higher than n — 1. Explicitly, the first few terms are given by 

F;f' {o) (X,t Q ;i,x i ) = (68) 
F A T> (1) (A,t ;t" &) = F^% =to>xi=0 (69) 
F^p t,(2) (A,t ;t,x J ) = F?^'^ \ t=tOjX i =0 + x*9ii^'^| t=t0jX i =0 + idoF^ t '^\ t=toX=0 (70) 



These equations again show explicitly the characteristic "mixing of orders" in the relation- 
ship between near-zone and far-zone fields. For example, equation ( 169]) demonstrates that 
the first-order near-zone electromagnetic field is a constant (in i, x l ) whose value is given 
by the value of the zerot/i-order far-zone field on the worldline at time to- Notice also 
that the zeroth-order near-zone external electromagnetic field vanishes, corresponding to 
the vanishing effects of the "external universe" in the near-zone limit. 

We turn now to the near-zone expansion of the electromagnetic self-field F^} f . The far- 

,">self 



zone expansion of was given above in eq. ( 1321) . Writing r = Ar and t = to + Ar in the 
first line of eq. (132|) . we obtain the corresponding near- zone series (dropping the remainder 



11 Note, however, that the far-zone multipole series was, in effect, defined with 7 taken to be at the origin 
of coordinates, whereas near-zone moments will be defined below relative to a A-dependent worldline, so 
the relationship between higher order far-zone and near-zone multipole moments will be complicated. 
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terms for convenience) 

N M 

F0\t, f, 6,<P) = J2Y1 ^ n ~ m - 2 (^u)nm(t0 + At, 9, <P) (71) 

n=0 m=0 
N M P 

= E E E A n+p ?r"- m - 2 (^) nmp (t ) (72) 

n=0 m=0 p=0 

where we have defined (J r flu )nm P = h {§i) P {^pLv)nm\t=t m the second line. Note that this 
expansion gives useful information at large f, in contrast to the far-zone expansion eq. (1321) . 
which gives useful information at small r. (Of course, these series contain precisely the same 
information, which is useful only at both small a = r = Xf and small (3 = X/r = 1/f.) Note 
that in the near-zone series, at nth order in A the highest possible combined positive powers 
of i and f is n, but arbitrarily high inverse powers of f can occur. Comparing with eq. ( |34l ). 



we can explicitly evaluate the first few orders of the near-zone series for F0 { at t = as 
follows 

fS^lw^ + O^) (73) 

Ff = O (^) (74) 

F!?- ( %. = I (-iferOn* - i*) + O (1) (75) 

F«"' (1) U = o(l) (76) 

^o lf,{2) |f=o = <? Q(%^) 2 nj + ^(ajn^ai + -a^rii + ~d rij + ^a^j + O f M (77) 

^7 lf ' (2) lf=o = -4^i + °(F) ' (78) 



where, in this equation and all following equations in this paper, we drop the bars on the 
coordinate components of barred quantities, i.e., it is understood that the components of 
any barred quantity is taken in barred coordinates. 
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B. Center of Mass and Body Parameters 



In order to determine the corrections to Lorentz force motion, we would like to define 
a notion of the "center of mass worldline" of the body at finite A. However, as discussed 
at the beginning of this section, this is a highly nontrivial notion for a charged body, since 
it is far from obvious how self-field contributions to the center of mass should be included. 
Fortunately, in order to define the first order correction to Lorentz force motion in the far 
zone, we only need to define a notion of center of mass to zeroth order in the near zone. As 
we shall now see, this can be done in a very simple and straightforward manner. However, 
if we wished to calculate higher order corrections, we would need to employ a much more 

n 

sophisticated notion of "center of mass" , such as given in [8J . 

To define the center of mass to first order in A, consider an arbitrary smooth, one- 
parameter family of worldlines 7(A) such that 7(0) is the worldline 7 of condition (i) of 
section II. Choose Fermi normal coordinates about 7(A), corresponding to working in the 
"rest frame" of the body. In these coordinates, the (flat) spacetime metric takes the form 
(see,e.g.,Q) 

g 00 = -1 - 2a 4 (A, t)x l - (a,(A, t)x 1 ) 2 + 0(r 3 ) 

g i0 = 0(r 3 ) (79) 

9ij "ij 7 

where a l (X,t) is the acceleration of 7(A) at time t. The components of the corresponding 
scaled metric are given by 

g 00 = -1 - 2\a i {to)x l - A 2 [2a i (t )x i i + (a^to)^) 2 + 25a i {t )x i ] + 0(A 3 ) 

9io = 0(A 3 ) (80) 

9ij "ij ) 

where we have defined a l (t) = a* (A = 0, t) and 5a l = <9a<2*(A = 0, t), and the overdot denotes 
the derivative with respect to t. Thus, a 1 corresponds to the (far-zone) acceleration of the 
unperturbed worldline 7(0), and 5a l corresponds to the perturbed (far- zone) acceleration of 
7(A) to first order in A. In accordance with the conventions stated at the end of the previous 
subsection, we have dropped the bars on the coordinate components of on the left side 
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of these equations. 

As before, we define the total body stress-energy tensor T a b by 




-iself 
ab 



(81) 



At zeroth order in the near-zone, the spacetime metric flHUj) is Minkowskian and the total 
body stress-energy is stationary (since all zeroth order near- zone quantities are stationary). 
Therefore, it is natural to define the total body mass at zeroth order in the near-zone by 



Since FjZ is a stationary, regular solution of Maxwell's equations that falls off at large 
f as 1/f 2 , the zeroth order self-field stress-energy falls off as 1/f 4 , so the integral defining 
m converges. Indeed, it is easily seen that m agrees with the total body mass at first 
order in the far- zone, eq. (|4"T|) . which we also denoted as m. Therefore, in particular, m is 
independent of t , a fact we will confirm in the next subsection by our near- zone analysis. 
We define the center of mass, A^ M , to zeroth order in the near zone by, 



The integrand in this equation falls off only as 1/f 3 , so the integral defining X l CM does 
not converge absolutely. However, the 1/f 4 part of Tqq (computed from equation (173|) ) is 
spherically symmetric, so the angular average of the 1/f 3 part of the integrand vanishes. 
Thus, the integral defining X % CM is well defined as a limit as R — > oo of the integral over a 
ball of radius R. 

The center of mass changes under a change of the origin of the near-zone coordinates 
x 1 . This corresponds to a first-order in A change of far- zone origin x l , or, equivalently, to 
a change in the choice of 7(A) to first-order in A in the far-zone. We therefore may define 
the perturbed motion to first-order in A in the far-zone by the demand that the zeroth order 
near-zone center of mass vanish, 



Since Xq M is defined relative to the Fermi (i.e., "rest frame") coordinate components of the 




(82) 




(83) 



^cm(^o) - . 



(84) 
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stress-energy, this demand corresponds to the requirement that an observer moving along 
the worldline assign his own position as the center of mass of the body at each time to- in 
the next subsection, we shall see by direct computation that this condition gives rise to an 
equation satisfied by the perturbed acceleration 5a l = d\a l \\ = o, thereby defining a correction 
to Lorentz force motion. 

We define the spin tensor to zeroth order in the near zone by 

S°S(t ) = -S j0 (t ) = J f^ 00 x j d 3 x (85) 

and 



S %3 (t ) = 2 J T m \xH s x. (86) 
Thus, when the center of mass condition is imposed, we have 

S 0j = . (87) 

We shall assume that the center of mass condition has been imposed in the following. As we 
shall see from eq. (11021) below, it follows from lowest order conservation of stress-energy that 
S* 2 - 7 is antisymmetric, S %:> = —S^ 1 . Note that the integral defining S 1 -* converges absolutely, 
since p^ 1 ^ ' falls off as 1/f 2 and _p^ clf '(°- ) falls off as 1/f 3 . The spin vector Si is defined by 

Si = -e ijk S jk . (88) 

To first order in the near-zone, the spacetime metric ( IHDl) is no longer Minkowskian, but 
it is stationary, with timelike Killing field d/dt. The stress-energy current T liv {d / dt) y is 
therefore conserved to first order, and it is natural to define the first order correction to the 
mass, 5m, by 



where the integral is independent of the spacelike hypersurface S. This yields 



(89) 

A=0 



5m(t ) = I T ( 1} _ d A x , (90) 



73^ 



t=Q 



where the center of mass condition has been used. Since the stress-energy decays only 
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as 1/f 3 on account of the "cross term" contribution of and , this integral 

does not converge absolutely. However, the explicit asymptotic form of T q (computed from 
equations (I73H76I) ) is 



+ , {91) 



mil) 

00 f=o " iTrr'"'" V ' Vr J 
from which it can be seen that the angular average of the 1/f 3 part vanishes, so that the 
integral defining 5m is well defined as a limit as R — ► oo of the integral over a ball of radius 
R. In view of eq. (1901) . it would be natural to interpret 5m as the first order correction to 
the rest mass of the body (including its electromagnetic self-energy). 
The total charge q(X) is given at finite A by 

q(X) = [ J a dZ aj (92) 

and is independent of the choice of hypersurface S. To zeroth order in the near-zone, the 
charge q is therefore given by 

q = J J^°d 3 x (93) 

It is not difficult to see that this agrees with the first order charge in the far-zone, which we 
previously also denoted as q. 

The zeroth order electromagnetic dipole moment tensor Q^ v of the body is defined by 

«-»>>-/ J"*'*. ^ 

together with Q j0 = -Q 0j . We shall see in the next subsection that the purely spatial 
components are automatically antisymmetric by lowest-order conservation of charge, so 
we have = Q [ ^ ] . The time -space components of Q^ v correspond to the electric dipole 
moment 

V 1 = Q° l (95) 
whereas the purely spatial components of Q^ v correspond to the magnetic dipole moment 

lh = -\e ijk Q jk ■ (96) 
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By eq.( |92|) . we see that the first order near- zone correction, 5q, to the charge is given by 

Sq = J JW°\t =0 d 3 x + a l Q 0l . (97) 

Note that the analog of the last term in this equation was absent from eq. (1901) on account 
of the imposition of the center of mass condition. Since q(X) is independent of to for all A, 
it follows immediately that 5q is independent of to- 



C. Derivation of Motion 



We now turn to the derivation of the perturbed equations of motion, which will be 
obtained directly from conservation of stress-energy, eq. (156|) . Writing T ab = + T^ lf as 
before, we obtain from eq. fl5T)|) 

V b f ah = J b F:f. (98) 

At zeroth order in the near- zone, the near-zone metric is flat (see eq. ( 180|) ). the external 
electromagnetic field vanishes (see eq. (1681)). and all zeroth order near-zone quantities are 
stationary (see eq. (1651)). Therefore, at zeroth order, we obtain 

&>f^ = 0, (99) 

where we again remind the reader that we have dropped bars on indices of barred quantities. 
By multiplying the time and space components of this equation by x l and integrating over 
space, we obtain 

J T^d 3 x = (100) 
J f^d 3 x = (101) 

Similarly, by multiplying the time and space components by x l x k and integrating over space, 
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we obtain 



f {0)l Q x k d 3 x = -J f (0)k Q x l d 3 x (102) 
f^x k d 3 x = 0. (103) 

Equation (11021) demonstrates the previously claimed anti-symmetry of the (lowest-order) 
spin of the body, defined by eq. (1861) . Equation fllQ3[) follows from an analogous anti- 
symmetry in % and k, combined with the symmetry of Tij in i and j. 
At zeroth order in the near-zone, conservation of charge-current yields 

Bt J (0)i = , (104) 

from which we can analogously derive the useful relations, 



J^d 3 x = (105) 
J^x j d 3 x = - [ J (0)J ¥ . (106) 



The second of these equations demonstrates the previously claimed anti-symmetry of the 
spatial components of the electromagnetic dipole tensor Q^ u defined by eq. (IMl) . 

At first-order in the near-zone, we see the first appearances of the acceleration (via 
the metric components flHOl ) and the external field (via equation fl69|) ). At first order, 
conservation of stress-energy, eq. (1981) . explicitly yields 

d°T^ + d%] ] + (jT$ - J^F*f = (107) 
9°f« + 8*1$ + a fc f ( 0) + a^f - J<f»»Fff = (108) 

where a« is evaluated at time t = t , and we have used eq. ( 1691 to replace F^'^ by 

pext — pext,(0) I 

r [iv — r p.v |i=i ,a; l =0- 

We shall now use these equations to derive evolution equations for the mass and spin to 
lowest order, as well as to determine the lowest order motion. For the evolution of mass, 
we integrate equation (11071) over space (with volume element d 3 x). Using the consistency 
condition ( 1661) . the first term evaluates to —dm(to)/dto. The second term vanishes by inte- 
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gration by parts, since T^' falls off like 1/r 4 . The third term vanishes by eq. (If 001) . and the 
fourth term vanishes by eq. fll05[) . Thus we obtain 

-^-m = , (f 09) 

in agreement with what we already found in our far-zone analysis given in section III. 

For the evolution of spin, we multiply eq. (11081) by x % and integrate over space. Again 
by the consistency condition the first term evaluates to — {l/2)dSik{to)/dtQ. The third term 
vanishes by the center of mass condition ( |84|) . and the fourth term vanishes by eq. ( 11031 ). 
However, the second and fifth terms do not vanish, yielding 



-— S lk = J T^d'x + CTiFff. (110) 
The antisymmetric part of this equation provides the desired evolution equation for the spin, 

^, = 2Q W ;. (in) 

The right side of this equation corresponds to the usual formula for the torque on an electro- 
magnetic dipole. It also should be noted that, since we work in Fermi normal coordinates, 
dSik/dt corresponds to the "Fermi derivative" of Sik, so eq. ( 11 111) automatically includes 
the "Thomas precession" (see eq. (H35p below). The symmetric part of (IllOp yields the useful 
relation 

'fV><Px = Cp {i Fffi. (112) 

To determine the lowest order motion, we integrate equation (11081) over space. The first 
term vanishes by the consistency condition and eq. (11001) . the second term vanishes since 



yield, 



T>j falls off as 1/r 3 , and the fourth term vanishes by eq. (I103p . The third and fifth terms 



mai = qF™ . (113) 



Thus, we reproduce the Lorentz force law previously derived by the far-zone analysis of 
section III. 

Finally, we derive an additional relation from first-order stress-energy conservation by 
multiplying eq. (I107p by x l and integrating over space. The first term vanishes by the 
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consistency condition (1661) and the center of mass condition fl84l) . The second term may be 
integrated by parts; the associated surface term vanishes because falls off as 1/f 4 . The 
third and fourth terms are simply expressed in terms of the spin (155]) and electromagnetic 
dipole (1§4"|) tensors. We obtain 

J f (1)j d 3 x = ^aiSV + Q ij F^ . (114) 

In a similar manner, at first-order in the near-zone, conservation of charge-current yields 

BojW + djJ^+af J™ = 0. (115) 

Multiplying by x % and integrating over space, we obtain 



/ 



jMtfx = -^Q 0i + aj Q ji . (116) 



We now shall derive an equation for the perturbed acceleration 5 a 1 (see eq. flHUl) ). as 
well as an evolution equation for the perturbed mass, 5m (see eq. (|90|) ). At second-order, 
conservation of stress-energy yields 

= a°f ( 2) + + a%] ] + SctT® + a%ft + 2a i x i B f^ + 2a i a; i f ( 0) - aWT^ 

- JW" (x%F°? + F^i + 5F°?) - jO-Kptf (117) 
= 5 f!£ + S*rff + a fc T « + a*T« + 5a k f^ + Sa'Tff + d k f^i+ itfffi 

- J®" (x%F«* + F^i+ SFtf) - J^F^, (118) 

where we have used eqs. (18*01) and (1701) and have written 5F^ = i 71 ^? |*=i ,a:*=o- To derive 
an equation for the perturbed acceleration 5a\ we evaluate (11 18[) at i = and integrate over 
space, applying the consistency condition (1661) to various terms, using the center of mass 
condition (18~4"1) . using eqs. (11 Oil) and (jl03j) . and using the definitions of the spin and dipole 
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tensors given in section HVBl We obtain 

= ^- J T^d 3 x + J fV'f^dfl + a k 5m + a 1 J T^d 3 x + 5a k m + j—Su + ~d% 
- Q^diF^f - 5F^ J J(°^rf 3 x - J JWPFffcPx . (119) 

The second term in this equation is the surface term arising from integration by parts of the 
second term in eq. (11181) . To evaluate this term, we compute the asymptotic form of T i 



(2) 



from equations ( 17311751) . obtaining 
1 q 2 



ij 



t=o Air r 2 



TiiUj [ 6(a k n k ) 2 + -a k a k + d ) + ( 4(a k n k ) 2 + -a k a k + -d + -a k n k 



A k 4. 1 



+ 0(1). (120) 



Since the integral over a sphere of the product of an odd number of normal vectors, n\ 

— (2) 

vanishes, only the terms in that are proportional to a 1 survive when computing the 
second term in eq. (11191) . We obtain 

J fVT^dJ] = 2 -q 2 ti (121) 

which can be recognized as the Abraham-Lorentz-Dirac self- force term (see eq. (Q above). 

Returning to eq. (11191) . we evaluate the first term using eq. (11141) . we evaluate the fourth 
term using eq. (II 12j) . and we use eqs. (11051) . (|97j) . and f II 1 6 j) to re-express the final two terms 
in terms of the electromagnetic dipole and the corrected charge. After some algebra, we 
obtain 

mSa, = -(5m)a l + (6q)F^ + q5F^ + p 2 a t - 1 -Q^d l F^ + ^ (o% - 2Q\,F^ ) , (122) 

where we have used Maxwell's equation d^F^ = to reexpress the fifth term on the right 
side. 

We also can derive an evolution equation for the perturbed mass, 8m, by integrating 
equation (11171) (evaluated at i = 0) over x l , and performing similar manipulations as above. 
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The result is, 

-±-5m = l -Q^d Q Fp - A ((Tiff ) • (123) 

Equation fll 22j) is our desired equation for the perturbed acceleration 6a\ As described at 
the beginning of the next section, this equation can be written as an evolution equation for 
the deviation vector, X a , describing the perturbation from Lorentz force motion. Appearing 
on the right side of eq. (11221) are the lowest order mass, m, the lowest order charge, q, 
the lowest order acceleration a\ the perturbed mass, 5m, the perturbed charge, Sq, the 
spin tensor, S^ v , and the electromagnetic dipole tensor, Q^ v . As previously derived, m, q, 
and Sq are conserved, and a 1 is given by the Lorentz force equation (11131) . The evolution 
of S^ v is given by eq. (II 111) , and the evolution of 5m is given by eq. (11231) . However, we 
cannot obtain an evolution equation for the electromagnetic dipole tensor Q^ v because the 
behavior of this quantity is "non-universal" . Indeed, it seems clear that by use of a "Maxwell 
demon" 12 , we could make Q^ u evolve with time any way that we wish. Therefore, additional 
conditions/assumptions beyond the Maxwell and matter equations CQ)-® and (J7J) would 
need to be adjoined to the equations we have derived in order to get deterministic evolution 
at this order. In appendix [B] we provide an example of a system with deterministic evolution 
by considering the case of a body with no electric dipole moment and a magnetic dipole 
moment proportional to spin. 

We can rewrite our equations of motion in covariant form as follows. In Fermi coordinates 
based on the worldline 7, the time direction (d/dt) a on 7 coincides with the unit tangent, 
u a to 7, so for example, can be replaced by M a F^ xt . Similarly, spatial components 
of quantities on 7 correspond to projections orthogonal to u a , so, for example, the spatial 
components correspond to the spacetime tensor (5 c b + w c M&)w a .F ( f* t = u a F^ ■ The 

spatial derivative, di, in Fermi coordinates can be expressed in terms of the spatial projection 
of the covariant derivative V M by using the following formula for the Christoffel symbol in 
Fermi coordinates, 

r Vl7 = aa u^u u - u a u^a u - u a a^u u . (124) 

Thus, diF^ 1 corresponds to (5 d c +u d u c )V ' dF^ +2a c u e F e ^u a ]. Finally, since, by construction, 
the Fermi coordinate vector fields {(d/dt) a , (d/dx l ) a } are Fermi transported along 7, it is 

12 Of course, the stress-energy of the Maxwell demon must be included in T™. 
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clear that taking the time derivative, "d/dto", of Fermi coordinate components along 7 
corresponds to taking the Fermi derivative of the corresponding tensors along 7. Here the 
Fermi derivative, Dp/dr, is given by the usual covariant derivative D/dr = w a V a along 7 
combined with a Lorentz boost in the plane defined by u a and a a . Specifically, for a scalar 
function / on 7 we have 

w f = (125) 

for a covector field f a on 7 we have 

^/« = ^/ fl + 2fa [A ]. (126) 
and for a tensor field f ab on 7 we have 

-jj: /«* = + 2 fa a [bU c ] + 2f c b a [a u c] . (127) 



In covariant notation, the evolution equations become 



m5a a = -a a 5m + 6[qF°?u b ] + a a u b u c Q hd F*? 

+ {9a + u a u b ) [p 2 ^ - ^Q cd V b F^ + ^ (a c S cb - 2u d Q c [b F$) } (128) 

S ab = 2 (<f a + u c u a ) (g d b + u d u b ) Q e [c F$ (129) 



In addition, we have u a S ab = (see eq. (187]) above). We remind the reader that these 
equations would have to be supplemented with a rule for the evolution of the electromagnetic 
dipole (arising from additional assumptions about the body) in order to comprise a "closed," 
deterministic system. 

By inspection, we see that eqs. ( 11281) and H130I) simplify if we redefine the perturbed mass 

by 

5m -f 5m = 5m- u b u c Q bd F°f . (131) 

This corresponds to adding the standard expression, — p ■ E ext , for electric dipole interaction 
energy (with E ey± the external electric field in the rest frame of the body) to the definition 
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of perturbed mass. Note, however, that we do not make any similar —jl ■ B ext adjustment 
to the definition of perturbed mass. With this re-definition, eqs. (11281) and (I130p become 



mda, 



a 



a a 5m+5[qFXu b } + {g a b + u a u b ) Sl q ^a b -^Q cd V b F^+^ {a c S cb - 2u d Q c [b F$) 

(132) 




We finally rewrite eqs. (1132}) . (11291) and (11331) in terms of the usual covariant derivative 
D/dr = u a V a along 7, obtaining 



where 5[ma a ] = m5a a + (5rh)a a . 

The first term on the right side of equation ( 11341) is simply the corrected Lorentz force. 
The first term in curly brackets is the Abraham-Lorentz-Dirac self force (see eq. (E]) above). 
The second term in curly brackets corresponds to the usual dipole forces familiar from elec- 
trostatics and magnetostatics. The final terms in curly brackets are not usually considered 
in elementary treatments of electromagnetic forces. We will show elsewhere [15J] that these 
terms are responsible for producing behavior in an orbiting charged body with spin and a 
magnetic dipole moment that is analogous to behavior found for spinning black holes in 
binary orbits in general relativity. The first term on the right side of eq. (11351) corresponds 
to the usual dipole torques familiar from electrostatics and magnetostatics. The second term 
is the Thomas precession. Finally, we note that the (modified) perturbed mass 8rh is not 
constant in time. However, in some cases — for example in the case discussed in appendix 
IB1 — there may be an alternative notion of mass that is conserved. 

Equations (1 1341) - ( 1X361) together with u a S ab = provide the complete description of the 
first-order deviation from Lorentz force motion, eq. (1481) . that can be derived from Maxwell's 
equations and conservation of stress-energy alone, and comprise the main result of this 
paper. 



5[ma a ] 



—5rh 
dr 



dr 
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V. SELF-CONSISTENT MOTION 



There are no mathematical difficulties that arise if one uses eqs fll34l) - fll36p to obtain the 
lowest order deviation from Lorentz force motion, eq. (TJHD- The lowest order quantities m, 
q, and are to be viewed as "given" , and the unperturbed worldline, 7, and its associated 
quantities u a , a a , and Da a / dr are to be viewed as having been determined by solving eq. (1481) . 
We wish to determine the perturbed worldline as well as the time evolution of the quantities 
S a b, Qab, 5q, and 5m. The perturbed worldline is described by a deviation vector, X a , on the 
background worldline 7, defined as follows: Consider the A-dependent family of curves, 7(A), 
introduced in subsection IIV Bl Parameterize each curve by proper time r. In arbitrary fixed 
(i.e., A-independent) coordinates x M , the family of curves is described by Z^{\\ r). Then the 
coordinate components of the deviation vector are given by X M = (<9Z M (A; t)/<9A)|a=o- It 
follows that 

5u a = ^-X a (137) 
dr 

and 

D D 2 

5a a = —5u a = — X a . (138) 
dr dr 2 

Furthermore, since SF^ = (J^ (F^iX) | 7 (a))) | a=0 we have 

5F:? = X c V c F^ + F a e b xt ' (1) , (139) 

where is the first order perturbation of the external field arising in far-zone pertur- 

bation theory, evaluated on 7. Again F^'^ should be viewed as "given" (and, indeed, it 
would normally be assumed to vanish). It can then be seen explicitly that eqs. (1 134f) - f TT36l) 
together with eq. f!137j) and D5q/dr = comprise a system of linear, first-order ordinary 
differential equations for (X a , 5u a , S a b, Q a b, 8q, 5m). If supplemented by an appropriate evo- 
lution equation for these equations have a unique solution for any given initial values of 
these quantities. The initial Sab also must satisfy the center of mass condition S a bU a = 0; 
this condition on S a & is then preserved by the evolution equation (I135p . The solutions to 
these equations are not plagued by the type of pathologies that occur for solutions to the 
Abraham-Lorentz-Dirac equation (Q. 

Nevertheless, the above equations do not give a fully satisfactory description of the devi- 
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ation from Lorentz force motion due to self- force, spin, and electromagnetic dipole moment. 
Even if all of the terms on the right side of eq. (11341) are small compared to the Lorentz 
force at all times, the cumulative effects of these terms over time should eventually make 
the deviation vector, X a , become large, at which point the description of motion as a linear 
perturbation of a single, fixed Lorentz force trajectory cannot be accurate. The solution to 
this difficulty is not to go to higher order in perturbation theory but to realize that although 
the deviation from a single Lorentz force trajectory may eventually become large at late 
times, the local deviation from some Lorentz force trajectory should be small at all times. 
As discussed in much more depth in [9[ , to implement the description of motion as a locally 
small deviation from some (varying) Lorentz force trajectory, we need to find a system of 
"self-consistent perturbative equations" that satisfies the following criteria: (1) They must 
have a well posed initial value formulation. (2) They must have the same number of de- 
grees of freedom as the first order perturbative system, so that a correspondence can be 
made between initial data for the self-consistent perturbative equation and the first order 
perturbative system. (3) For corresponding initial data, the solutions to the self-consistent 
perturbative equation should be close to the corresponding solutions of the first order per- 
turbative system over the time interval for which the first order perturbative description 
should be accurate. 

An obvious first attempt to find such a self-consistent perturbative equation would be to 
simply "delete the 5's" from eqs (1134jl - (1136p . i.e., take the corrected equations of motion to 
be 

qF^u b + {g a b + u a u b ) ^<f^a b - \o^h^t + £ {a c S cb - 2u d Q c [b F$) } 

(140) 

2 (g a c + u a u c ) (g b d + u b u d ) Q e [c F$ - 2a c S c[a u b] (141) 
lQ*§*T + 2Q a b F b Ta [c u^ (142) 

together with u a S ab = 0. In a sense, this is the simplest modification we can make to the 
Lorentz force equation (HHl) that takes into account the first-order perturbative effects of 
eqs. (I134p - (11361) . However, it is easily seen that although eqs. (I140p - (11421) are ODEs and 
thus have an initial value formulation, these equations violate criteria (2) and (3) above. 



ma„ 



, ^>ab 
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D . 



dr 



-m 
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Specifically, the role of the acceleration, a a } changes from that of a zeroth order "background" 
quantity in eq. fjl34j) to that of an unknown dynamical variable in eq. (11401) . Consequently, 
the terms involving Da a / dr on the right side of eq. (11401) make this equation be of higher 
differential order than the corresponding perturbative equation (11341) . Thus, more initial 
data is required for eqs (I140j) - (jl42p than for eqs (I134l) - (ll36p . in violation of criterion (2). 
Criterion (3) also is violated. Indeed, when Q ab , and S a b vanish, eq. (11401) reduces to the ALD 
equation Qj, so even in this simple case, eqs (I140p - (11421) have solutions that differ drastically 
from eqs (I134p - (ll36p . It is also worth noting that even when q and Q ab vanish (so there 
are no electromagnetic effects), if S a b is nonvanishing, eq. (11401) is still of higher differential 
order than eq. (I134p on account of the D(a a S a b)/dr term. The resulting unphysical degrees 
of freedom of the system eqs (I140p - (ll42p together with u a S a b = give rise to solutions with 
"helical motions" in addition to the expected solutions with inertial motion. 

However, this difficulty can be overcome by the reduction of order procedure described 
in the Introduction. We replace Da a /dr in the two terms in which it appears on the right 
side of eq. (I140p by the right side of eq. ( flTT) . and we replace a a on the right sides of 
eqs. (1 1401) -( fbiiZl) by the right side of eq. ( flOl) . The resulting system of equations — which we 
will not write out explicitly here — still corresponds to Lorentz force motion corrected by 
the perturbative effects of eqs. (I134p - (I136I) . but the resulting equations now have highest 
derivative terms of exactly the same form as eqs. (I134p - (ll36p . Consequently, criterion (2) 
holds, and we believe that criterion (3) also holds. Thus, we believe that the reduced order 
form of eqs (ll40l) -f fT4"2l) supplemented by u a S a b = comprise a fully satisfactory system of 
self-consistent perturbative equations that provide a description of motion that takes into 
account the lowest order effects of self-force, spin, and electromagnetic dipole moments. For 
comparison with what can be found in textbooks on electromagnet ism, we note that in the 
non-relativistic limit 13 and in ordinary vector notation, eq. (I140p takes the form 

F = q(E + vx B)+ ^q 2 ^ + Pi S7 E l + fi&B 1 

+ ^- (s x a — ax E — px B] , 
dt V / 

13 To obtain this equation, we drop all terms in eq. (| 140|) that are quadratic or higher order in velocity, v, 
as well as terms linear in v that are multiplied by q 2 , p, fl or S. 
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(143) 



where we have written F = ma and it should again be understood that the terms in a and 
da/dt on the right side of this equation should be eliminated by reduction of order. 

Finally, as discussed in [9] , there is no reason to expect self-consistent perturbative equa- 
tions satisfying criteria (l)-(3) to be unique. Although the reduced order form of eqs (I140p - 
(11421) supplemented by u a S a b = appears to be a fully satisfactory system of self-consistent 
perturbative equations corresponding to the corrections to Lorentz force given by eqs (EEMj)- 
), other choices of self-consistent perturbative equations are possible. In particular, the 



following system of equations arises naturally from the analysis of the motion of extended 
bodies [8(: Define the "force" f a and "torque" n ab = n\ a u in terms of the 4-momentum P a 
and and electromagnetic dipole Q ab by 



f a = 2 -{q/Mfh a b P d {P c V d F^ - qftFgFf) - ^Q bc V a F b e c xt (144) 



n ab = 2Q c [a F^ , (145) 



where 



and 



M = V-P a P a , (146) 



hob = 9ab + PaPb/M 2 . (147) 

(The first term in eq. (11441) can be recognized as a reduced order form of the Abraham- 
Lorentz-Dirac force.) Then the tangent, f y a , to the center-of-mass worldline (normalized so 
that P a 7 a = —M) is given in terms of P a , Q ab and the spin tensor S ab by 



- P n P b lM S ^c t {P c -n c d P d /M) + Mf b ] 
M la -P a - n ab P /M M2 _ • 1 1 l b > 



The evolution equations for P a and S ab are 

P a = qFlTi h + fa (149) 

S ab = 2P [a 7 b ] + n ab . (150) 

These equations are to be supplemented by P a S ab = 0, which, if imposed as an initial 
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condition, can be shown to be preserved under evolution (see [16j). At the level of what 
is known purely from lowest order perturbation theory as derived in section IV, there is 
no reason to prefer the system (I144j) - (|15U|) to the system (I140j) -( TT^^T) . However, in cases 
where the self-force term in eq. (I144p is negligible, the system (I144p - (11501) yields an exact 
conservation law when there is a Killing field that Lie derives — a feature that would 
hold exactly for test body motion. Therefore, it appears that in this respect the system 
ffT44l)- frT50|) is superior to the system ffM -fTM. 

VI. SUMMARY 

We have given a rigorous and systematic treatment of particle motion in classical electro- 
magnetism. We considered a one-parameter-family of solutions to the Maxwell and matter 
equations ([I])-© and ([7]) containing a body that "shrinks down" to zero size, mass, and 
charge according to the scaling assumptions of section [Til We found that the lowest-order 
description of the body is that of a point particle (1291) and (T46l) moving according to the 
Lorentz force law (1481) . The first-order corrections to this motion — including self-force, dipole 
force, and spin force effects — are then given by our rigorous perturbative result (I134I) - (I136I) . 
Finally, we addressed the issue of finding a self-consistent perturbative equation associated 
with our perturbative result. We argued that the naive self-consistent perturbative equations 
(ll40l) -( fT42l) can be modified by "reduction of order" to provide appropriate self-consistent 
perturbative equations. In the case of negligible spin and electromagnetic dipole moment, 
this reduces to the reduced-order ALD equation. 

Acknowledgements This research was supported in part by NSF grants PHY04-56619 and 
PHY08-54807 to the University of Chicago. 

APPENDIX A: SCALING PROPERTIES OF THE SELF-FIELD 

This appendix derives the two properties of F*^(\, t, x % ) that were claimed to hold at the 
end of section [TTJ for families of current distributions satisfying eq. ( fl9l) . 

We first show that the retarded field F^ lf (A, t, x % ) satisfies the scaling relation, eq. (TSUI) . In 
global inertial coordinates (t, x l ), the Lorenz-gauge vector potential for the retarded solution 
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with source J M is given by 

Af(\,t,x>) 



dV 



J„(X,t 












- X H 





(Al) 



where, in our case, we have ^(A,^,^ 4 ) = A _2 J M (A,t, [x l — z l (t)]/X) with J M smooth in all 
variables and of compact support in the spatial variables. Now introduce scaled coordinates 
x % = [x l — z l (t)]/\ and x H = [x H — z % (t)]/\. Using these definitions together with eq. ([19]) 
and the coordinate transformation y l = x H — x i , it is easily seen that 



A s ; [i (\,t,x*) = J d'y 

Since z l (t) is smooth, there must exist a smooth function V 1 such that 

z \t) - z%t - X\y j \) = X\y j \V%X\y j \). 



J^\t-X\f 


, x j + y j + [z j (t) - z j (t - A 


2/ fc |)]/A) 









(A2) 



(A3) 



For later reference, we note that since x % = z l (t) defines a timelike curve, it follows from 
the mean value theorem that \V l \ < 1. Substituting in eq. (1A2|) . we see that the integrand 
is smooth in (A, t, x l ) and that the convergence properties of the |/-integral are sufficient to 
allow us to interchange differentiation with respect to (A, t, x l ) with integration. It follows 
that in global inertial coordinates, A s ° lf smooth in these variables, i.e., ^4^ elf is of the form 



Af(\,t,x*) = A,(\,t, [xW(t)]/A). 



(A4) 



where ^(A, t, x l ) is smooth in all of its arguments. Differentiating with respect to the space- 
time variables, we find that F^(\,t,x l ) takes the form ( 1201) in global inertial coordinates. 
By the same type of argument as given in sect. [TTl it then follows that F^ lf (A, t, x l ) takes 
the general form (|20|) in arbitrary coordinates for which V M t is timelike, as we desired to 
show. 

We now derive a certain falloff behavior for F^jr outside of the body, specifically that, in 
terms of the parameters 

a = \x l - z\t)\, P = X/a, (A5) 
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we have 



XF^f = (3 2 F, l/ (t,a,(3,n i ) (A6) 



at least for a and /3 in a sufficiently small neighborhood of (0,0), where is smooth in 
all of its arguments and 

rj = [x* - z\t)]/a , (A7) 

so that \n l \ = 1. An analogous statement applied to one-parameter families of metrics 
was used in |9j in the analysis of gravitational self-force; as explained further in 9], this 
corresponds to what is needed to justify the use of matched asymptotic expansions. Matter 



distributions (if any) were not explicitly considered in [9], so the analogous behavior of 
metrics in that reference was an assumption. Here, we prove that the electromagnetic 
analog follows from a very natural form for the source, namely eq. ( fl9l) . 

We start, again, with eq. flAll) for the Lorentz-gauge vector potential A s ^ lf (A, t, x l ) in 
global inertial coordinates. Initially, we assume that a > 0, so that (3 > and n % are well 
defined. We define 

w * = x ' 1 -z\t- \x j -x' ] \) (A8) 



and we use eq. ( 1A3I) to write 

z\t - \x j - x' j \) = z\t) - \x j - x'^V^t, \x j - x' j \) . (A9) 



We also write 



w i = \w { = aPw' 1 . (A10) 



The integrand in eq. (IA1I) is proportional to J^(\,t — \x l — x H \, tip). Since J M is of compact 
support in the last argument, there exists a positive constant D such that only points 
satisfying \w l \ < D contribute to the integral. Consequently, we may assume that \w l \ < D 
in the following. 

Equations (lA8l) -( TAT0l) together with the above definitions of a and n l yield 



x l - x 11 = a{n l - (3w t ) + - x' 3 \V l {t, \x k - x' k \) . (All) 
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Thus, taking the squared magnitude of both sides of this equation, we obtain 

y 2 = \ a {n i -Pw i )+yV i {t } y)\ 2 . (A12) 

where we have written 

y = 1^ -x^l . (A13) 

Treating eq. (1A12|) as though it were a quadratic equation in y (even though V 1 depends on 
y), one finds that 

(\V k \cos6 + y/l - \V k \ 2 + \V k \ 2 cos 2 , (A14) 



a\ 


n % - 




1 - 




2 



where cos 6 is defined by 



\v>\ m e = v '} n ! - (Ai5) 



and V 1 is evaluated at (t, y). Since y appears on both the left and right sides, eq. OA 140 is 
not a solution for y but rather is a relation that must be satisfied by the variables y, t, a, 

(3, n\ and w\ 

As already noted below eq. (1A3I) . we have \V l \ < 1. Furthermore, since \w t \ < D, if we 



restrict to f3 < 1/D, then \n l — f3w l \ > 0. Consequently, for (3 < 1/D, the right side of 
eq. (1A14I) is of the form of a times a smooth, nonvanishing function of (t,?/,n l , (3w l ). We 
may therefore solve this equation for a, obtaining a solution of the form 

a = yH(t,y,n i ,pw i ) (A16) 



where H is a smooth, nonvanishing function of all of its arguments. 

Thus far, we have assumed that a > 0, which implies y > given our restrictions that 
|u> l | < D and /3 < 1/D. However, eq. (1A141) makes perfectly good mathematical sense — and 
H in eq. (1A16|) remains smooth — if we extend the domain of y and a to include 0. (Indeed, 
eq. (1A14I) makes sense — and H remains smooth — even for (unphysical) negative values of y, 
(3, and a.) We may therefore apply the inverse function theorem at y — a — to conclude 
that in a sufficiently small neighborhood of a — 0, y can be written as 

y = aL(t, a,n\(3w j ) , (A17) 
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where L is a smooth, positive function of its arguments. 
Rearranging (IA8I) . we have 



x n = a(3w l + z l (t- aL(t,a,n\(3w k )) . (A18) 

We now return to eq. flAlj) and make the coordinate transformation x' % — ► to 1 . Using 
eq. (1A18|) . we see that the Jacobian of this coordinate transformation takes the form 



dx' 1 



{aP) 3 W(t,a,p,n\w j ) , (A19) 



dutf 

where W is a smooth function of its arguments. Substitution into (lAlj) finally shows that 



The denominator here has been shown to be strictly positive over the range of w l for which 
the numerator is non-vanishing. Both L and W are smooth and the integral is being taken 
over a compact region. The integral is therefore smooth in t, a, (3, and n l , i.e., there exists 
a smooth A^(t, a, j3, n l ) satisfying 

Af = PA, (A21) 
for sufficiently small a and j3. Differentiating to obtain the field strength shows that 

F^f = (P/a)^ u (t,a,P,n l ) (A22) 

for some smooth in all of its arguments. Since A = a/3, it follows immediately that 
AF^ lf satisfies eq. ( 1A6I) . We have derived eq. (1A6|) in global inertial coordinates, but it is 
readily seen that this form is preserved under an arbitrary smooth transformation to new 
coordinates for which V^t is timelike. 



APPENDIX B: "ELECTRON" MOTION 

It has been stressed in subsection IIVCI that the perturbed equations of motion that can 
be derived using only Maxwell's equations and conservation of total stress-energy do not 
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constrain the time evolution of the electromagnetic dipole moment Q ab . Thus, in order to 
have a deterministic system, the equations we derived in subsection IIVCI would have to be 
supplemented by an evolution law for Q ab arising from the particular matter model under 
consideration. A particularly simple and interesting case is to assume that for some (time 
independent) constant C we have 

Qab = —CS a b. (Bl) 

From eqs. (|14|) and (|15|) together with our center of mass condition u a S ab = 0, it is easily 
seen that this is equivalent to having a vanishing electric dipole moment of the body and a 
magnetic moment proportional to its spin, i.e., 

p a = 0, fi a = CS a . (B2) 

Since we have derived an evolution equation for S a b, eq. (IB1I) serves to fix the time evolution 
of Q a b as well. As discussed at the beginning of section [V] eqs (I134|) - (ll36p together with 
eq. (IBip and u a S ab = comprise a well-posed deterministic system for determining the 
perturbed motion. Note that since the electric dipole moment vanishes, we have 5rh = 5m 
(see eq. f TTail ). 

Since many elementary particles satisfy eq. ( IBip . we may view the self-consistent pertur- 
bative system eqs. (I140p - (11421) together with eq. (IBip as providing a simple, classical model 
for the motion of an elementary particle, such as an electron. To write these equations in 
a more familiar form, it is useful to decompose into rest-frame electric and magnetic 
fields via 

F a c b xt = 2u [a E b] + e abcd u c B d . (B3) 
The evolution equation (11411) for the angular momentum can then be written as 

^S a = -e abcd u b n c B d + 2u [a S b] a b . (B4) 

Since fi a = CS a , the right side is orthogonal to S a , so 

^(S a S a ) = 0, (B5) 

i.e., the spin vector does not change its magnitude under time evolution. Clearly, /i a also 
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does not change its magnitude, i.e., the body has a "permanent" magnetic dipole moment. 
The evolution equation (11421) for mass reduces to 



where we have dropped the "hat" on m since there is no distinction between rh and m when 
the electric dipole moment vanishes. Although m does not remain constant, it is clear that 
we may define a new quantity 



that is conserved. Note that the standard expression for interaction energy, —\i a B a , is being 
subtracted from m in this equation — rather than added to m — to define m*. 

Note that the failure of the "rest mass" m to be constant resolves a paradox concerning 
what one is taught in elementary physics courses: On one hand, one is (correctly) taught 
that an external magnetic field can "do no work" on a body, so a body moving in an external 
magnetic field cannot gain energy. On the other hand, one is (also correctly) taught that a 
magnetic dipole released in a non-uniform external magnetic field will gain kinetic energy. 
Where does this kinetic energy come from? Equation (1B6I) shows that it comes from the 
rest mass of the body. 1 

We now consider the motion of the body, as described by eq. (1140)) . Using eq. (IBID , we 
obtain 



where the terms involving acceleration on the right side are understood to be eliminated 
by reduction of order. 15 If q ^ 0, it is conventional to write C in terms of a dimensionless 

14 By contrast, for a "permanent electric dipole", i.e., Q a i, — 2it[ a p h i with DpQab/dT = 0, it follows from 
ea. (|142p that the original mass m = rh + p a E a is conserved. The kinetic energy gained by an electric 
dipole released in a non-uniform electric field comes from the work done on the body by the electric field. 

15 In particular, a term proportional to eb c dfd c n d E^ has been set to zero since a a is to be replaced by qE a /m. 




(B6) 



m* = m + fi a B a . 



(B7) 




(B8) 
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g-f actor defined by C = qg/2m, in which case 



CE a -a a = 



(9 - 2)q 



E 



•a 



(B9) 



2m 



where we have again used a a = qE a /m for terms on the right side. In the case of an electron 
or muon, g is very nearly 2, so the last two terms on the right side of eq. flB8j) should be 
negligible. However, g can be very different from 2 for composite particles, in which case 
the last two terms in eq. ( 1B8I) need not be negligible. 

Another point worth noting with regard to both eq. (I140p and eq. ( 1B8I) is that the 
Lorentz-Dirac self-force is not the only term involving Da^/dr. A "jerk" force also occurs 
in conjunction with the spin of the body in the second line of eq. (IB8I) . (Note that this 
force is orthogonal to the ALD force.) For a body of charge q equal to the charge of the 
electron e, and with spin of order h, the ratio of the magnitude of the "spin jerk" term to 
the ALD term is of order (g — 2)h/e 2 = (g — 2)/a. Thus, if g ^ 2, the "spin jerk" term can 
dominate over the ALD term. However, the effects of the "spin jerk" term probably do not 
accumulate significantly over time. 

It also interesting to ask when the ALD self-force becomes comparable to the Lorentz 
force. In this situation, it is likely that higher order corrections to the self-force will be 
large, and our perturbative equations are unlikely to be reliable. It is easily seen that for an 
external electromagnetic field that oscillates with frequency to, the condition that the ALD 
self-force is much smaller than the Lorentz force is simply 



It seems unlikely that a classical description of electron motion would be possible in any 
case if this condition were violated. 

Another criterion for the validity of our perturbative equations is that the change in mass 
under time evolution — which is a perturbative effect — be small compared with the mass, i.e., 
that 



q 2 uj <C m . 



(BIO) 



For an electron, this reduces to 



(Bll) 



n\\B\ <C m 



(B12) 
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If g ~ 1 and \q\ = e, this reduces to 

\ B \ «^~(4x 10 13 G) ( -J^-] 2 (B13) 
1 1 eh V ; V500keV/ V ; 

Magnetic fields expected to exist near some neutron stars are expected to violate this bound. 
It seems unlikely that any presently known classical equations of motion can adequately 
describe the behavior of electrons when this bound is violated. 
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